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Abstract 

The branching ratio of the lightest Higgs decay into bb, rf and cc is sensitive 
to supersymmetric effects. We include in this work the effects of CP phases on the 
Higgs decays. Specifically we compute the deviation of the CP phase dependent 
branching ratio from the standard model result. The analysis includes the full one 
loop corrections of fermion masses including CP phases involving the gluino, the 
chargino and the neutralino exchanges. The analysis shows that the supersymmet- 
ric effects with CP phases can change the branching ratios by as much as 100% for 
the lightest Higgs decay into bb and rf with similar results holding for the heavier 
Higgs decays. A detailed analysis is also given of the effects of CP phases on the 
Higgs decays into cc. The deviations of Rb/ T an d Rb/ C from the standard model 
result are investigated as a possible signature of supersymmetry and CP effects. 
Thus a measurement of the decays of the higgs into bb, rf and cc may provide 
important clues regarding the existence of supersymmetry and CP phases. 



: Permanent address of T.I. 



1 Introduction 



Is is known that supersymmetric contributions can significantly affect the Higgs 
decays into bb, rf and cc[l]. We compute here the effects of CP phases on these 
arising from soft breaking parameters. Thus, for example, the simplest supergrav- 
ity unified model mSUGRA[2], whose soft breaking sector is defined by the param- 
eters too, mi, A , tan/5 (where m is the universal scalar mass, mi is the universal 
gaugino mass, A is the universal trilinear coupling and tan/3 =< H 2 > / < Hi > 
where H 2 gives mass to the up quark and Hi gives mass to the down quark and the 
lepton) can accommodate two CP violating phases which can be chosen to be the 
phase of the Higgs mixing parameter \i and the phase of A . Extended versions of 
mSUGRA including nonuniversalities can accommodate more phases. In this anal- 
ysis we will consider the more general case of the supersymmetric standard model 
(MSSM) which allows for several CP phases. Thus, for example, we will allow the 
gaugino masses to be in general nonuniversal so that rhi = \rhi\e % ^(i = 1,2,3). 
In the analysis we use the results recently obtained regarding the effects of CP 
phases on the third generation quark and lepton masses [3]. The analysis of Ref.[3] 
extends the analyses of Refs. [4, 5, 6] where large corrections to third generation 
masses and specifically to the b quark mass, but without inclusion of CP phases, 
were found. The corrections to the quark and lepton masses are of considerable 
importance in the analysis of Yukawa unification [7, 8, 9]. In this paper we utilize 
these corrections to study their effects on the Higgs decays into bb and rf. 

The analyses which include CP phases must be constrained by the experimental 
upper limits on the electric dipole moment for the electron and for the neutron 
which are very stringent. Thus for the electron the current experimental limit for 
the magnitude of electron edm is d e < 4.3 x ICT 27 ecm[fO] while for the neutron 
it is d n < 6.5 x 10~ 26 ecm[ll]. There is a similar stringent limit on the edm 
of the Hg" atom, i.e., dn g < 9 x I CT 28 ecm[I2]. A variety of techniques have 
been discussed in the literature to achieve consistency with experiment [13, 14, 15, 
16, 17, 18, 19, 20, 21]. Specifically one finds that it is possible to accommodate 
large CP phases and still have consistency with the edm constraints either via the 
cancellation mechanism[16, 18, 20, 21] or via large phases in the third generation 
sector [19]. We note in passing the analyses of Refs. [20, 21] also include H^ 99 edm 
constraint and show that this constraint along with the electron and the neutron 
edm constraint can be satisfied in the presence of large phases. Of course, if the 
phases are large they will affect a variety of low energy phenomena. Such effects 
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have been investigated on a variety of processes. These include effects on the Higgs 
sector [22, 23], on — 2 [24], on collider physics [25, 26], in B physics and in flavor 
violation [27, 21, 28, 29] and a variety of other low energy phenomena. The number 
of phenomena investigated is rather larger and a more complete list can be found 
in Ref.[30]. 

The outline of the rest of the paper is as follows: In Sec. 2 we give a brief 
description of the basic formalism needed for the evaluation of the full one loop 
effects on the decays of the Higgs into quarks and leptons including the effects of CP 
phases. In Sec. 3 we give the effective low energy interaction of the b quark with the 
lightest Higgs including loop corrections. The branching ratios BR(H 2 — > bb) and 
BR(H 2 — > tt) are also computed. In Sec. 4 we take the limit of vanishing phases 
and compare our results to previous analyses. In Sec. 5 we give a numerical analysis 
of the deviations of the branching ratios from the standard model predictions 
because of supersymmetric effects and discuss the sensitivity of these deviations 
to CP phases. Conclusions are given in Sec. 6. In Appendix A we extend the 
results of Sec. 4 to include H 2 — > cc decay and compute the supersymmetry and 
CP effects on AR b / c . Extension of the results to decays of the Higgs bosons Hi 
and H 3 is given in Appendix B. 

2 The basic formalism 

At the tree level the b quark couples to the neutral component of Hi Higgs boson 
while the coupling to the H 2 higgs bosons is absent where 

(*.) = ($). = (3) w 

Loop corrections produce a shift in the H® couplings and generate a non vanishing 
effective coupling with H®. Thus the effective Lagrangian would be written as 

-C eff = (h b + 5h b )b R b L H° + Ah b b R b L H°* + H.c. (2) 

where the star on H®* is necessary in order to have a gauge invariant £ e //- The 
same analysis holds for the tau-lepton sector where h T , 5h T and Ah T are used in 
the Lagrangian involving trTl. 

The quantities Shf and Ahf receive SUSY QCD and SUSY electroweak con- 
tributions. They are calculated in Ref . [3] on mass corrections to lepton and quark 
masses. In the analysis carried out there one finds that the couplings are generally 
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complex due to CP phases in the soft breaking terms. Electroweak symmetry is 
broken spontaneously by giving vacuum expectation value to H® and H®. Thus 
one finds 

Inserting in and H% one finds 

-C m = M b b R b L + H.c. (4) 

where 

M b = — =-{1 + — + —— tan/5} (5) 

Here M& is complex because and Ah b are both complex. We carry out a 
redefinition of the b quark field 

b = e*™l/, tan Xfe = ^^ (6) 

ReM b 

After the redefintion of Eq.(6) the mass term reads 

-C m = m b b' R b' L + H.c. (7) 

where m b is real and positive and b' is the physical field and m b and h b are related 
by 

h b = y^[(l + Re 6 -^ + Re^t a nPf + (Im^ + Im^ tan/?) 2 ] \ (8) 
vi h b h b h b h b 

The above can be approximated by 

V2m b 1 Sh b Ah b 

h b ~ — — — , A 5 = Re— + Re-— tan/3 9 

vi (1 + A b ) h b h b 

From now on we will drop the prime on b and we will assume that we are already 
in the basis where the b field is the physical field for the b quark. The other terms 
in Eq. (2) which have (f>i,tpi,(/>2 and ip 2 dependence will produce the interaction 
between the b-quark and the mass eigen states of the Higgs fields H^. Next we 
introduce the basis <f>i, fa, i^id, i>2D where 

ipin = sin flip! + cos Pip 2 , ^2D = - cos fiipi + sin fifa (10) 
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In this basis the field ip2D is the would be Goldstone field. By considering one 
loop contributions to the Higgs boson masses and mixings from top-stop, bottom- 
sbottom, chargino and neutralino exchanges with CP violating phases, the mass 
eigen states of the Higgs fields would be mixed states of CP even and CP odd 
states. Thus the mass eigen states Hi (i=l,2,3) are related to the 4>i,4> 2 ,ip\D a s 
follows [22, 23, 31, 32] 







( 01 ^ 




= R 


02 


\H 3 J 







In the absence of CP phases there is no mixing between the CP even and the CP 
odd Higgs and the Higgs mass 2 matrix consists of a 2 x 2 matrix for the CP even 
fields and a one element for the CP odd Higgs. Solving for the eigen values in this 
case one may choose the first eigen value to be the heavier mass and the second 
eigen value to be the lighter mass. In previous analyses where we considered the 
effects of CP phases on the Higgs sector [23], we used the convention that in the 
limit of vanishing CP phases one has H 1 — > H,H 2 — > h, H 3 — > A [23]. In the 
present analysis we continue to use the same convention. Thus the lightest Higgs 
boson field corresponds to H 2 in our notation. 

3 The Decays H 2 — ► bb and H 2 — ► rf including 
effects of CP phases 

In this section we study the decay of H 2 Higgs into bb and rf . Consider the b 
quark first. The effective interaction of the b quark with the Higgs mass eigen 
states H 2 is given by 

-C\ nt = b[C s b + t l5 C[]bH 2 

Cf = Cf cosx 6 -Cf sin X& 
Cf = Cf sin Xb + Cf cos Xb 
v^Cf = Re(h b + 5h b )R 21 + [-Im(h b + 8h b ) sin (3 
+Im(Ah b ) cos(3}R 23 + Re(Ah b )R 22 
V2C b = -Im(h b + 5h b )R 21 + [-Re(h b + 5h b ) sin (3 

+Re(Ah b ) cos (3]R 23 - Im(Ah b )R 22 (12) 

In deriving C\ nt above, we redefined the quark field b such that h b is real and is 
given by Eq.(8) while the quark mass M b is complex before field redefinition. In 
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Figure 1: One loop contribution to the bottom quark mass involving exchange of 
gluino, charginos and neutralinos in the loop. 

Ref. [31] the authors redefined the fields such that h b is complex and is given by 
their Eq.(4.5) while the quark mass is real and positive. To compare now our 
Eq.(12) with Eq.(4.10) of Ref. [31] and Eq. (77) of Ref. [6], we have the phase 
Xb of the quark mass apparent while h b in Cf and is real. In Eqs.(4.11, 4.12) 
of Ref. [31] and Eq.(78, 79) of Ref. [6] the phase Xb is absent because of their 
choice of phases which makes hb + Sh b + Ahb tan (5 a real and positive quantity and 
our Eq.(6) would lead then to a vanishing Xb- However h b in these expressions is 
complex and this will compensate for the absence of Xb- So both methods, the one 
given here and that of Refs.[31] and [6], lead to the same phenomenolgical results. 



The r lepton has similar interactions. Thus 

-Cl nt = flC? + t l5 C?}bH 2 (13) 

where and are given by Eq. (12) with the transposition b — > r. We are 
interested in the ratio of the branching ratios for H 2 — > bb and H 2 — * fr, i.e., 

_ BR(H 2 -> 66) 
Rb/T ~ BR(H 2 - fr) (14) 

Using the interactions given in Eqs. (12) and (13) we can estimate this to be 

= 3^^(l + a;) (15) 
rib Db J-t 



Here 



T/M M^ Mr+( Mh_^ mf (16) 

1VI H 2 1V1 H 2 
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where 



A s f = A s f cos X f - A p f sin Xf , A P = A s f sin X j + A P cos X f 
A b f = [(1 + R e ^l)R 21 - Im—±R 23 sin/3 + Re—±R 22 + Im—±R 23 cos/3] 



1 ' ^A^-lJ^R^P + Re^Rv + Im^ 
h f h f h f hf 

A P = [-I m 5 -hlR 2l - (1 + Re 6 -^)R 23 sin/3 - Im^R 22 + Re^R 23 cos/3](17) 
h f hf h f h f 



Sf = [(1 + Re 6 -^ + Re^ tan/3) 2 + (Im^- + Im^ tan/3) 2 ] (18) 
h f h f h f h f 

and (1 + a;) is the QCD enhancement factor and is given by[33] 



rv rv 

A + uj) = 1 + 5.67— + 29.14^ (19) 
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so that + ~ 1.25 for a s ^ 0.12. CP phases enter in Eq. (15) through Tf and 
Sf since these depend on the phases through the couplings Ahf, 5hf and through 
the matrix element Rn. Now the CP dependence of Ahf, and 8hb is significantly 
different from the CP dependence of Ah T and 5h T on phases since, for example, 
corrections to the b quark mass involve gluino contributions and are thus sensitive 
to the phase £3 while the corrections to the r lepton mass do not depend on this 
phase. Thus one can expect a sensitive dependence of the ratio Rb/ T on the phases. 
We may compare the above result to the result from the Standard Model. Here 
one has 

(R b /r)sM = 3(4)( m !"! m b ( 1 + <*) (2°) 

By identifying rrih with m# 2 we define the shift due to supersymmetric effects 
including the effects due to CP phases as follows 

AR b/T = Rb/T - {Rb/r)sM (21) 

Rb/T 

As pointed out in Ref.fl] this ratio can be used to distinguish the SM Higgs from 
the lightest SUSY Higgs. The same holds here except that we also take into 
account the effects of CP phases. We will discuss the effects of the CP phases 
on ARb/ T numerically in Sec. 5. It is also of interest to analyze the ratio of the 
branching ratios H 2 — > bb and H 2 — > cc. Thus define 

BR(H 2 ^bb) 
Rb/C ~ BR(H 2 -+ cc) (22) 
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By repeating the analysis as for the previous case we get 

ml „S'T h 



R b/C = ^ 2 t^p^ (23) 
mi S b Tl 



- c 



where 



where 



A s c = if cos X c - if sin xc, A p c = if sin X c + if cos X c 
if = [(1 + rJMr 22 - Im^R 23 cos (3 + Re^R 21 + Im^R 23 sin/5] 

if = [- Im 5 -hlR 22 _ (i + Re 5 -^)R 23 cos (3 - Im^R 21 + Re^R 23 sin 0\ (25) 

^ C ' 'C c c 

# = [(i + # e ^£ + i? e ^ cot/5) 2 + (Jm^ + Jm^ cot/5) 2 ] (26) 

c c c 

We note that the prime is used on T' c and S' c for the c quark case since it has 
different couplings. The Ah c and Sh c can be deduced similar to the analysis of 
Aht and Sht given in Ref. [3]. Explicit results are given in Appendix A. An analysis 
similar to the above but for the decay of the heavier Higgs bosons Hi and H 3 are 
given in Appendix B. 

4 Limit of Vanishing Phases 

We compare our results now with previous analysis of Ref. [1] by taking the limit 
of vanishing phases. In this limit we set 

Xf = 0, Im(8h f ) = = Im(Ahf) 
Rn = cos a, R\ 2 = sin a 
R 2 i = — sin a, R 22 = cos a 

Rl3 =R23 = (27) 

Neglecting the masses of b, c and r relative to the masses of the Higgs bosons, one 
obtains from Eqs. (15)- (18) the result 

_ o m 2 b (1 + 4-e b /tana) 2 (1 + e' T + e T tan/5) 2 
Kb/T ~ m? (l + e' b + e fe tan/5) 2 (1 + e' T - e T / tana) 2 [L + u) [ ^ } 
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where e b)T and e' b r are defined so that 



e b:T = Ah b:T /h b:T e' bT = 5h b>T /h b>T (29) 

On neglecting the r correction and e' b correction this result agrees with Eq.(ll) of 
the first paper of Ref. [1]. Similarly in the limit of vanishing phases one finds for 
R b /c the result 

o m l u . m 2 i 1 + 4, ~ e b cot af (1 + e' c + e c cot (3f 

Rb/c = — ^(tanatanp) t^tz ; ^ (30) 

' m 2 c y ' ' (l + e' b + e 6 tan/?) 2 (l + e^-e c tana) 2 v ; 

Again neglecting the e' b and e' c terms the above result agrees with Eq.(19) of the first 
paper of Ref. [1] . We note, however, that in general the e' term is not necessarily 
negligible. 



5 Numerical Analysis 

In this section we discuss the size of the supersymmetric corrections including CP 
phases to AR b / T and AR b / c discussed in Sec. 3. We carry out the analysis in MSSM. 
Since the general parameter space of MSSM is rather large we shall limit our selves 
to a more constrained set for the purpose of this numerical study. We shall use for 
our parameter space the set m^, m , m 1 / 2 , \Aq\, tan/3, 6> M , cxa , £i, £2 and £3. The 
parameter \i is determined via radiative breaking of the electroweak symmetry. 
The other sparticle masses are obtained from this set using the renormalization 
group equations evolving the GUT parameters from the GUT scale down to the 
electroweak scale. We then use the sparticle spectrum generated by the above 
method to compute the supersymmetric corrections to the Higgs decays. As men- 
tioned above the phases chosen here for the numerical analysis are £ 1; £ 2 , £3, (Xao 
and Op. This choice provides a convenient set of phases because of their appear- 
ance in the expressions for Rij, Ahf and Shf. By examining these expressions one 
finds linear combinations of the phases listed above that enter the analysis. Such 
linear combinations have been classified for EDMs in Ref. [17]. Thus not all the 
phases are independent and absorbing one of these phases by redefinition of the 
other phases in these linear combinations will not change the numerical analysis, 
i.e, the effect of CP phases on Higgs decay. However, for the numerical analysis 
here it is more convenient to work with the set listed above. 

In Fig. 2 we give a plot of AR b / T for the decay of Higgs H 2 as given by Eq. (21) as 
a function of the phase £3 for values of tan (3 ranging from 5-50. The analysis shows 



8 



a very sharp dependence of ARb/ T on £3. Thus £3 affects not only the magnitude of 
ARt)/ T but also its sign depending on the value of £3. A similar analysis holds for 
AR^ T where R^L is given by Eq.(54). Fig. 3 gives a plot of AR^ T as a function 
of £3 for all the same inputs as for Fig. 2. In this case also one finds a very sharp 
dependence of AR^} on £ 3 . The analysis of AR^? follows a similar pattern as 
AR^ T and is not exhibited. As discussed in Sec. 4 the quantity ARb/ c is also of 
interest where Rb/ C is given by Eq.(30). In Fig. 4 we give a plot of ARb/ c as a 
function of £3 for all the same inputs as in the analysis of Fig. 2. Here also we find 
a very substantial dependence of ARb/ c on £3. Next we analyse the variations of 
ARb/ T and Ai4/ C as a function of 9^. In Fig. 5 a plot of ARb/ T is given as function 
of for various values of £3. Similar plots for AR^} and for ARb/ c are given in 
Fig. 6 and Fig. 7. Together Figs. 5, 6 and Fig. 7 show a rapid dependence of 
ARb/ T , AR^} and of ARb/ c on 9^ in addition to their sharp dependence on £3. 

The analysis given in Fig. 2, Fig. 3 and Fig. 4, and also in Figs. 5, 6 and Fig. 7 
include CP effects arising from both supersymmetric QCD and supersymmetric 
electroweak effects. Next we study the dependence of the branching ratios on the 
electroweak CP phase £ 2 - In Fig. 8 we give an analysis of ARb/ T for the decay of 
Higgs H 2 as a function of £ 2 for values of m , mi ranging from 200 GeV -400 GeV. 
Here we see that the dependence of AR b / T on £ 2 is not as strong as it is on £3. 
Nonetheless the £2 dependence of ARb/ T is still quite substantial as the variations 
in AR b / T can be as much as 40%. In Fig. 9 we exhibit the dependence of AR^ T 
on £2 and one finds that once again the dependence though not as strong as the 
£ 3 dependence is still quite substantial and one can get variations of as much as 
40% over the full range of f 2 . The analysis of AR^ T is very similar and is not 
displayed. In Fig. 10 the dependence of AR b / c on £ 2 for the H 2 decay is given. 
Again although the dependence of AR b / c on £ 2 is not as strong as on £ 3 it is still 
quite substantial as one finds that AR b / c can vary about 20% over the full range 

of 6- 

In Fig. 11 we give a plot of ARb/ T as a function of tan/3 for three different 
inputs which for the largest value of tan (5 satisfy the edm constraints including 
the Hg" edm constraint. (The H^ 99 edm constraint translates into C Hg = \d% — 
d% - 0.012df J < 3.0 x 10- 26 cm[20, 21].) The lower set of curves are for the cases 
when phases are included while similar upper curves are without phases. First one 
finds the effect of phases is so large as to not only affect the magnitude of AR b / T 
but also affect its sign. Further, one finds that AR b / r is strongly dependent on 
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tan (3 for cases with and without phases. An identical analysis for AR b } is given 
in Fig. 12 where AR^ r is plotted as a function of tan (3 for exactly the same set 
of inputs as in Fig. 11. One finds once again a very similar behavior in that both 
the magnitude and the sign of Ai?^_ are affected and also one finds a very strong 
dependence on tan f3. The analysis of Ai?^_ is very similar to that of Ai?^ and 
is not given. In Fig. 13 we give a plot of A/4/c as a function of tan/5 again for the 
same identical inputs as in Fig. 11. In this case also the effects of phases though 
smaller than in the case of ARb/ T are still substantial in that the CP phase effects 
can be up to 30%. Further, one finds about 20% variation due to variations in 
tan p. 

Table 1: EDM constraints for the points of tan/3 = 50 of Figs. 11 , 12 and 13. 



Case 


\d e \e.cm 


\d n \e.cm 


C Hg cm 


(i) 


1.67 x 10- 2Y 


1.59 x 10~ 2Y 


1.18 x 10~ 2Y 


(ii) 


6.05 x 10~ 28 


3.47 x 10~ 27 


1.29 x 10~ 26 


(iii) 


2.14 x 10~ 27 


8.90 x 10~ 28 


1.25 x 10" 26 



6 Conclusions 

In this paper we have computed the effects of CP phases on the Higgs boson decays. 
In supersymmetry after spontaneous breaking one is left with three neutral Higgs 
bosons which in the absence of CP phases consist of two CP even Higgs and one 
CP odd Higgs. With the inclusion of CP phases the Higgs mass eigen states are 
no longer CP eigen states but rather admixtures of CP even and CP odd states 
when loop corrections to the Higgs boson masses are included. Further, inclusion 
of loop corrections to the quarks and lepton masses are in general dependent 
on CP phases and these effects can be very significant for the case of the b quark 
mass. Additionally inclusion of CP dependent loop effects on the quark and lepton 
masses modify the vertices involving the quarks and leptons and the Higgs mass 
eigen states and these modifications also affect the Higgs boson decays. In this 
paper we have computed the effects of CP phases on the decays of the light and the 
heavy Higgs boson decays to bb and to rf. Specifically we computed the deviation 
of i4/r from the standard model prediction in the presence of CP phases. We find 
that these effects can be rather large. Thus ARb/ T can vary by as much as 100% 
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or more due to the variation in the 9 ^ and £3 while variations with £2 are relatively 
small although still substantial. A similar analysis was also carried out for Ai?5/ C 
and one finds that the variations of ARb/ c with phases are also substantial. The 
analysis presented here points to the possible detection of CP phases via accurate 
measurement of the branching ratios provided we have sufficient constraints on the 
remaining SUSY input parameters from other experiment. Thus the decays of the 
Higgs if measured with sufficient accuracy may provide a signal for the presence 
of both supersymmetry and CP phases. 
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Appendix A: Analysis of A.h c and Sh c . 

Following the same technique as given in Ref..[3] we give below explicit expressions 
for Ah c and 8h c . Thus the c quark mass at the Z scale is given by 

m c (M z ) = h c (M z )^= sin/3(l + A c ) (31) 
where A c gives the loop correction to the c quark mass m c . and A c is given by 

A c = (Re^cot(3 + Re 5 ^) (32) 
An analysis of Ah c to one loop order gives 

A^ = -EE ^e^m § F^D^D c2j fiml, m|, m\) 
i=ij=i 671 

222 m , 

- E E E 9 2 H* clJ {U* kl D: it - k s U; 2 D: 2 Mk c V: 2 D sl3 )^J(m 2 xt ,ml,ml) 

i=l j=l k=l 

222 m +m + 

-T,T,E9 2 K*Au; 1 D: lk -kM; 2 D: 2k }(k c v; 2 D^ 
t=i j=i k=i iD7r 1 3 

224 mo 
- E E E lF* cij {a ck D cl3 - l ck D c2j }{(3l k D* cli + a ck D* c2i }-^f(m 2 x m|, m|) 

i=ij=ik=i iD7r 
4 4 2 m Y om Y o 
- E E E lA\ 3 {a C3 D clk - l cj D c2k }{(3* ci D* clk + a a D: 2k }^-^f(ml,m 2 x0 ,m 2 x o)(33) 



, / , 167T 2 
=1 j=l k=l 



where 



g 2 m c X 4k 

OLck 



2mw sin j3 



(3 ck = eQ c X'* k + -i— X 2 * k (T 3c - Q c sin 2 6 W ) 

COSVw 

7 c fc = eQ c X lk — — X 2k (34) 

COS t>w 
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and where Q c — \ and T 3c = |. In the above D ci j is the matrix that diagonalizes 
the c squark mass 2 matrix and q are the c squark mass eigen states so that 

2 2 

cl = J2 D ciiCi, c R = ^2 D C 2iCi (36) 
1=1 1=1 

Similarly D si j is the matrix that diagonalizes the s squark mass 2 matrix and 
are the s squark mass eigen states so that 

2 2 

s L = D siiSi, sr = Y1 D s2i Si (37) 
i=i i=i 

Further, U and V are the matrices that diagonalize the chargino mass matrix and 
the matrix X diagonalizes the neutralino mass matrix and the elements of the X' 
matrix are defined by 

X' lk = X ik cos 9 W + X 2 k sin Q w , X' 2k = —X lk sin 9 W + X 2k cos 9 W (38) 

where we are using a notation of Ref. [3]. F ci j and H ci j are defined by 

■"' <J UZ {(I ~ Urn 2 9 w )D* cli D clj + \ sin 2 9 w D* c2i D c2j } cos/3 



^2 2cos^ LV 2 3 wy c1 ' cy ' 3 



and 



H dj _ ^ u 1 , 1,2^ ^* n 1,2. 



~ 2cos^ {( ~2 + 3 Sin e w)D* ali D slj - - sin 2 9 w D* s2i D s2j } cos/3 



[D* sU D slj + D* s2i D s2j ] - » ' " ^ (40) 



2M w cos/3 L sil ilJ s ^ AZJJ 2M^cos/3 
while and are defined by 



2^' vl~ 2^" (41) 



where 



= \l 2 Ui2Vjl 



9Q"n = \ K(gx*2j - gfx& + (< — i)] (42) 
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Finally the function / (m 2 , m 2 , m 2 ) is given by 



f(m 2 ,m 2 ,m 2 ) = - — - — ^ — ^ ^(m 2 m 2 ln^ + rr^m^ln^ + m 2 m 2 /n^§)(43) 



(m 2 — m 2 )(m 2 — m 2 )(m 2 — m 2 ) 3 m 2 1 mf 1 3 m 2 ' 



for the case i ^ j and 



f(m 2 ,m 2 ,m 2 ) = ( ^ ? _\ _ ( m ^ w ^ + ( m 2 - m 2 )) (44) 
for the case i=j. Similarly for 5h c we find the result 

Sh c = - E E m c 2 8 , m|) 

i=Jj=l J7r 

2 2 2 j }"! -f 

222 m +m + 

"EEE JC&U^D^ - k s U; 2 D: 2k }(k^D slk )^-^f(ml, m\ p m 2 +) 

i=l j=l k=l 

2 2 4 mo 
"EEE 2S CJ - i {o! cfc D c i : ,- - lckD c2j }{P* ck D* cU + a cfcJ D* 2 J— ^f(m 2 , ,, m? , m|) 
i=i j=i fc=i i07r 

4 4 2 m ofiij 

-EEE sr^-K^ - 7c A^H/mv + «^: 2 J-|-^/K, m % <°X 45 ) 
i=i j=i fc=i 

where E ci j and G ci j are given by 

7i = " f^M^cy + ^sin 2 ^: 2 ^ c2 ,}sin/5 

^ Ki^cy + D: 2i D c2] ] - ^^D: 2i D c1j (46) 



2M w sin/T M J CM 31 2M w smj3 



and 



G '" (J1[Z {(4 + I sin 2 9w)ir Ai D.u - \ sin 2 OwD^D^} sin/3 



^2 2cos0 w lv 2 ' 3 3 



while Cij and are defined by 



l 9m s y 



and 



T2 = -dr^ s "- Q '^- R ^ < 48 > 
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where 



R'ij = 2M^^ x *i x *^ + ^x* 2i x* 2j - pi*(x* 3i x: 3 + X* 4i X* 3j )} (50) 



Appendix B: Decays of Hi and H 3 Higgs bosons. 

We can repeat the same analysis for the heavy Higgs bosons Hi and H3. The 
effective interaction of the b quark with the Higgs mass eigen state Hi is given by 

-C b mt = b[Cf + i^C^bHi 
Cf = Cf cos xt - C b sin Xb 
C b = C b smxb + C b cos Xb 
V2Cf = Re(h b + 5h b )Ru + [-Im(h b + 5h b ) sin (3 
+Im(Ah b ) cos /3}R 13 + Re(Ah b )R 12 
V2C h p ' = -Im(h b + 5h b )Ru + [-Re(h b + 5h b ) sin (3 

+Re{Ah b ) cos l3]R n - Im(Ah b )R 12 (51) 
The r lepton has similar interactions. Thus 

-a mt = f [Cf + h,C p t '\tHi (52) 

and Cf and Cf can be obtained from Eq. (51) by the interchange b — > r. For 
the Hi decays one finds 



Hl _ BRjHi -> 6ft) 
Using the interactions given by Eqs. (51)and (52) we dertmine this to be 



(53) 



where 



where 



R "^W {1+U) (54) 



rf . = ( ^M)i (4? + ( ^4 )i(4? (55) 



Af = Af cos X f - Af sin X/ , = Af sin X/ + if cos X f 

7 4' r/ ^ <^/\ ^ t „ „ Ahf „ A/if „ 

= 1 + Re—^)Rn - Im—±Ri 3 sin /? + Re—^R U + Im—±R U cos /3 
1 h f h f h f h f 

Af = [-Im^-Rn - (1 + Re^-)Ri 3 sin (3 - Im^-R 12 + Re^^Ri 3 cos /?](56) 
h f h f h f h f 
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The interaction that governs the H 3 decay is 

-C\ nt = b[C s b " + t l5 C b p "]bH 3 
c b" = Cf cos Xb ~ C b p " sin Xb 
C b =C b smxb + C b cosxb 
V2C b " = Re(h b + 5h b )R 31 + [-Im(h b + 5h b ) sin (3 
+Im(Ah b ) cos (3]R 33 + Re(Ah b )R 32 
V2C b " = -Im(h b + 5h b )R 31 + [-Re(h b + 5h b ) sin (3 

+Re(Ah b ) cos (3}R 33 - Im(Ah b )R 32 (57) 

Similarly the r lepton interaction with H 3 is given by 

-C7 mt = f\Cf + h,C p "}tH 3 (58) 

and Cf and C p " can be obtained from Eq. (57) by the interchange b — > r. For 
the H 3 decays one finds 



>H3 BR(H 3 -» bb) 
hb ' T BR(H 3 -> fr) 

Using the interactions given in Eqs. (57) and (58) we dertmine this to be 

?H 3 _ o m b T b A 
h l T mlS h T?s 

where 



R »/r ~ BR(H 3 - fr) (59) 



<t = 3^|^£(l + -) (60) 



where 



= + (%^K"> 2 (61) 



Af = Af cosxf- Af sinx/, Af = Af sin^/ + Af cos^/ 

7<?" r/ ^ $hf .^ 5hr „ „ „ Ahf ^ Ahf „ 

Af = 1 + Re—±)R 31 - Im—±R 33 sm(3 + Re—±R 32 + Im—^R 33 cos/? 



Af = [-Im^R 31 - (1 + Re^)R 33 sin (3 - Im^R 32 + Re^R 33 cos/5](62) 

Alf /If /if Alj 
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Figure 2: Plot of AR b / T for the decay of H 2 Higgs as a function of the phase £ 3 . 
The other input parameters are: = 200 GeV, mo = mi = 200 GeV, & = .5, 

2 

£ 2 = -7, 9^ = .1, olao — 1-0, and \A \ = 4. The curves in descending order at 
the point £ 3 = tt correspond to tan/3 = 5,10,20,30,50. All angles here and in 
succeeding figure captions are in radians. 
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Figure 3: Plot of AR b / T for the decay of H 1 Higgs as a function of the phase £ 3 . 
The other input parameters are: = 200 GeV, mo = mi = 200 GeV, £i = .5, 

2 

£ 2 = -7, 6^ = .1, (Xaq = 1.0, and \A Q \ = 4. The curves in descending order at the 
point £3 = 7r correspond to tan (3 = 5, 10, 20, 50, 30. 
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Figure 4: Plot of AR b / c for the decay of H 2 Higgs as a function of the phase £ 3 . 
The other input parameters are: = 200 GeV, mo = mi = 200 GeV, £i = .5, 

2 

£ 2 = -7, 6^ = .1, chaq = 1.0, and \A Q \ = 4. The curves in descending order at the 
point £3 = 7r correspond to tan (3 = 5, 10, 20, 30, 50. 
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Figure 5: Plot of AR b / T for the decay of H 2 Higgs as a function of the phase 
The other input parameters are: = 200 GeV, m = mi = 300, £i = .5, f 2 = 

2 

dA = 1.0, \A \ = 5, and tan (3 = 20. The curves in descending order at 6^ - 
correspond to £3 = 0, .75, 1.5, 2.25, 3 
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Figure 6: Plot of AR b / T for the decay of Hi Higgs as a function of the phase 
The other input parameters are: = 200 GeV, m = mi = 300, £i = .5, f 2 = 

2 

dA = 1.0, \A \ = 5, and tan/? = 20. The curves in descending order at 6^ - 
correspond to £3 = 0, .75, 1.5, 2.25, 3 



24 




e^(radian) 



Figure 7: Plot of AR b / c for the decay of H 2 Higgs as a function of the phase 
The other input parameters are: = 200 GeV, m = mi = 300, £i = .5, f 2 = 

2 

dA = 1.0, |v4 | = 5, and tan (3 = 20. The curves in descending order at 6^ = 
correspond to £3 = 0, .75, 1.5, 2.25, 3 
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Figure 8: Plot of ARb/ T for the decay of H 2 Higgs as a function of the phase 
£2- The other input parameters are: uia = 200 GeV, £1 = .5, £3 = 0, 0^ — .1, 
cxa = 1.0, \Aq\ = 5, and tan/3 = 50. The curves in ascending order correspond to 
m = mi= 200, 250, 300, 350, 400 GeV. 
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Figure 9: Plot of ARt/ T for the decay of H\ Higgs as a function of the phase 
£2- The other input parameters are: uia = 200 GeV, £1 = .5, £3 = 0, 0^ — .1, 
cxa = 1.0, \Aq\ = 5, and tan/3 = 50. The curves in ascending order correspond to 
m = mi= 200, 250, 300, 350, 400 GeV. 
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Figure 10: Plot of ARb/ c for the decay of H 2 Higgs as a function of the phase 
£2- The other input parameters are: uia = 200 GeV, £1 = .5, £3 = 0, M = .1, 
oia = 1-0, \Aq\ = 5, and tan/3 = 50. The curves in ascending order correspond to 
m = mi = 200, 250, 300, 350, 400 GeV. 
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Figure 11: Plot ARb/ T for the decay of H2 Higgs as a function of tan/3. The 
three lower curves are for the cases (i), (ii) and (iii) as follows: (i)m A = 200 GeV, 
m = m i = 300 GeV, A = 4, a Ao = 1, & = .5, f 2 = -659, £ 3 = .633, 9^ = 2.5 
(solid) ; (ii) = 200 GeV, m = mi = 555 GeV, A = 4, a Ao = 2, £1 = .6, 
£ 2 = .653, £ 3 = .672,0,, = 2.5 (dot-dashed); (iii) m A = 200 GeV, m Q = mi = 480 
GeV, A Q = 3, a Ao = .8, ^ = .4, £ 2 = .668, ^ 3 = .6, 9^ = 2.5 (long-dashed). The 
edm constraints including the Hjj 39 are satisfied for the above curves at tan /? = 50 
as shown in Table 1. The three similar curves in the upper half plane are for the 
three cases above when the phases are all set to zero. 
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Figure 12: Plot ARb/ T for the decay of H\ Higgs as a function of tan/3. The 
three lower curves are for the cases (i), (ii) and (iii) as follows: (i)mA = 200 GeV, 
m = m i = 300 GeV, A = 4, a Ao = 1, & = .5, f 2 = -659, £ 3 = .633, 9^ = 2.5 
(solid); (ii) = 200 GeV, m = mi = 555 GeV, A = 4, = 2, £i = .6, 
6 = -653, £ 3 = .672,^ = 2.5 (dot-dashed); (iii) m A = 200 GeV, m = mi = 480 
GeV, A = 3, a Ao = .8, & = -4, 6 = -668, £ 3 = .6, M = 2.5 (long-dashed). The 
edm constraints including the Hjj 39 are satisfied for the above curves at tan /? = 50 
as shown in Table 1. The three similar curves in the upper half plane are for the 
three cases above when the phases are all set to zero. 
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Figure 13: Plot ARb/ c for the decay of H2 Higgs as a function of tan/3. The 
three lower curves are for the cases (i), (ii) and (iii) as follows: (i)mA = 200 GeV, 
m = m i = 300 GeV, A = 4, a Ao = 1, & = .5, f 2 = -659, £ 3 = .633, ^ = 2.5 
(solid); (ii) = 200 GeV, m = mi = 555 GeV, A = 4, = 2, £x = .6, 
6 = -653, £ 3 = .672,0,, = 2.5 (dot-dashed); (iii) m A = 200 GeV, m = mi = 480 
GeV, A = 3, a Ao = .8, £1 = .4, £ 2 = -668, £ 3 = .6, # M = 2.5 (long-dashed). The 
edm constraints including the Hjj 39 are satisfied for the above curves at tan /? = 50 
as shown in Table 1. The three similar curves in the upper half plane are for the 
three cases above when the phases are all set to zero. 
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